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l/"~) ' Ultracold quantum gases offer a wonderful playground for quantum many- 

^^ , body physics [1], as experimental systems are widely controllable, both stati- 

cally and dynamically. One such system is the one-dimensional (ID) Bose gas 
on a ring. In this system binary contact interactions between the constituent 
bosonic atoms, usually alkali metals, can be controlled in both sign and mag- 
k> , nitude; a recent experiment has tuned interactions over seven orders of mag- 

5^ ■ nitude [2], using an atom-molecule resonance called a Feshbach resonance. 

Thus one can directly realize the Lieb-Liniger Hamiltonian (LLH) [3, 4], from 
the weakly- to the strongly-interacting regime. At the same time there are 
a number of experiments utilizing ring traps [5]. The ring geometry affords 
us the opportunity to study topological properties of this system as well; one 
of the main properties of a superfluid is the quantized circulation in which 
the average angular momentum per particle, L/N, is quantized under rota- 
tion [6]. Thus we focus on a tunable ID Bose system for which the main 
control parameters are interaction and rotation. We will show that there is 
a critical boundary in the interaction-rotation control-parameter plane over 
which the topological properties of the system change. This is the basis of our 
concept of metastable quantum phase transitions (QPTs). Moreover, we will 



2 Understanding Quantum Phase Transitions 

show that the finite domain of the ring is necessary for the QPT to occur at 
all because the zero-point kinetic pressure can induce QPTs, i.e., the system 
must be finite; we thus seek to generalize the concept of QPTs to inherently 
finite, mesoscopic or nanoscopic systems.* 

Specifically, we will show that past the critical boundary the phase of a 
Bose-Einstein condensate (BEC) can wind and unwind continuously and can 
therefore no longer be considered a superfluid in the usual sense, despite a 
single mode still being macroscopically occupied [7] . Since the interaction and 
rotation will be in units related to the ring size, another way to look at this 
critical boundary is that the superfluidity in this example is a mesoscopic effect 
and might vanish in the thermodynamic limit. Using mean-field theory in 
the form of the celebrated Gross-Pitaevskii or nonlinear Schrodinger equation 
(NLSE) , the presence of the critical boundary can be tied to the appearance 
of new stationary solutions, dark solitons, which bifurcate from a uniform 
superflow or plane-wave branch. Dark solitons, which take the form of a 
density notch in BECs, can have zero density, in which case they are called 
black and have a node, or a region of decreased density without a node, in 
which case they are called gray. Gray solitons have a characteristic phase 
structure which allows for the topological winding and unwinding, as we will 
show. As the density contrast of a gray soliton tends to zero, i.e., as it 
becomes shallower, its phase structure, density, and energy approach those 
of uniform superflow. Arrays of regularly spaced dark solitons are called 
dark-soliton trains. Dark solitons, black and gray, singly and in trains, have 
been observed in BEC experiments, beginning with [8, 9] and, more recently, 
in [10, 11]. Although dark-soliton trains are robustly stable in quasi-lD mean- 
field theory,^ in practice they are found to decay in experiments. This has 
been attributed first to use of 3D geometries, in which dark soliton nodal 
planes decay via the snake instability into vortex-anti- vortex pairs; in this case 
mean-field theory provides an excellent description of the ensuing dynamics. 
However, even in quasi-lD experiments dark solitons are found to decay; the 
main cause is thought to be thermal [13] and/or quantum fluctuations [14, 15]. 

The study of quantum fluctuations, and the use of dark-soliton decay as a 
smoking-gun signal for fluctuations, has a long track record. It turns out that 
dark solitons can also be observed in optical fibers with anomalous dispersion 
and a Kerr nonlincarity [16]. The governing equation is again the NLSE, and 
dark solitons decay due to periodically spaced fiber amplifiers. This effect is 
called Gordon-Haus jitter [17], and can be formulated as a ID random walk 
with the accompanying diffusion equation. In BECs, thermal effects can play 
a similar role [18, 13]. However, at very low temperatures and for a quasi-lD 
mean-field one is left with only quantum fluctuations as the cause of dark soli- 



*See also Chap. 27 for a discussion of finite-size effects on QPTs in the context of nuclei. 
tQuasi-lD means that the mean field is confined but the underlying binary scattering 
problem remains 3D. See [12] for details. 
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ton decay. Thus a thorough investigation of such fluctuations was undertaken 
by Dziamarga et at, starting with [14]. Using the Bogoliubov-de Gennes 
equations (BdGE) and modified version thereof, they found that although 
dark solitons delocalized, due to the anomalous or Goldstone mode, they did 
not in fact decay. That is to say that averaging over many measurements 
one finds the soliton fills in, but for a single measurement a dark soliton is 
found somewhere in the condensate. However, the BdGE are limited to weak 
interactions. Recently one of us studied this problem in the more strongly in- 
teracting regime, placing the soliton in an optical lattice to enhance quantum 
fluctuations and using matrix-product-state (MPS) algorithms to follow the 
dynamics."' 1 We found that a dark soliton does indeed decay at zero temper- 
ature, based on three main pieces of evidence: density-density correlations, 
inelasticity in soliton-soliton collisions, and a close comparison to Bogoliubov 
theory via an analysis of modes of the single-particle density matrix. It has 
recently been demonstrated that it is possible in the weakly interacting regime 
for the filling-in of density-density correlations to be consistent with dereal- 
ization rather than decay [19, 20]; however, the strongly-interacting regime 
remains an open question. Thermal fluctuations in dark solitons have also 
been used to study the dynamics of phase transitions in the Kibblc-Zurck 
mechanism [21]; see Chap. 3. 

We will show that dark-soliton trains on a ring are in fact the weakly- 
interacting limit of yrast states [22, 23], and that it is really the yrast states 
that play the key role in our finite-size metastable QPT when analyzed in 
the full many-body problem beyond mean-field theory. Yrast, a Swedish term 
originally used in nuclear physics which can be translated as 'dizziest,' refers 
to the lowest energy state for a given angular momentum. We will identify 
Lieb's Type II or hole excitations with the yrast states, solving the long- 
standing question of the physical interpretation of this solution branch of the 
LLH [24]. The uniform superflow states which the soliton trains bifurcate 
from will be replaced with the more general center- of-mass rotation (CMR) 
states, a special class of yrast states for which L/N is an integer. In the 
NLS mean-field theory, the phase winding number comes out of a calculation 
of the circulation; out of the mean-field regime, such terminology becomes 
questionable if not meaningless. Thus the whole concept of uniform superflow 
and of superfluidity breaks down; yet our QPT persists as an accurate picture 
of an abrupt transition in even the strongly-interacting system's properties. 
In fact, we can take this description all the way to the Tonks- Girardeau (TG) 
limit in which one obtains an exact solution to the problem via the TG Bose- 
Fermi mapping [25]. At the end of our chapter, we will come back to the 
connection between yrast states and quantum fluctuations of dark solitons. 

Finally, we wish to make a brief comment on dimensionality. Although 
there is formally no supcrfluid phase transition in ID, nevertheless for an 



' Specifically, wc used time-evolving block dccimation(TEBD); sec Chap. 23. 



4 Understanding Quantum Phase Transitions 

interacting finite system there is in practice still a sharp crossover to a macro- 
scopically occupied mode at a critical temperature. One should keep in mind 
that bounded and/or finite domain reduced-dimensional systems often have 
distinctly different properties from their infinite counterparts, for instance 
graphene, which is unstable on the infinite domain due to thermal fluctua- 
tions, but, as we know now, not only stable on finite domain but has many 
important technical applications [26]. Likewise, this is an exciting time for 
the study of ID systems, not so much for the possibility of revolutionizing 
computing as with graphene, but for revolutionizing our understanding of 
quantum many-body physics. ID systems have some special advantages for 
that purpose. First, there are established exact solution methods such as the 
Bethe ansatz for special but nevertheless important systems, including the 
focus of the present study, the LLH. We will take advantage of the plethora 
of ID theoretical techniques in our study in order to bring to bear multiple 
lines of evidence on our problem. Second, experiments in ultracold quantum 
gases [27] are allowing for the thorough investigation and testing of funda- 
mental ID questions, including Kolmogorov-Moser- Arnold (KAM) theory for 
intcgrable quantum systems, quantum quenches, and a host of other zero 
and finite-temperature effects. Third, there are exciting new developments in 
MPS numerical methods that allow one to follow entangled quantum many- 
body dynamics of ID systems, integrable or nonintegrable, as discussed in 
Chaps. 22 and 23. 

This chapter draws strongly from three papers [28, 29, 30], the foundations 
for which were laid by [31]; however, our previous work on dark solitons [32, 
33, 34, 15] also plays an important role in our discussion and our thinking. 
Our set of theoretical techniques for this study included the NLSE, the BdGE, 
exact diagonalization in a truncated angular momentum basis, the finite-size 
Bethe ansatz equations, and the TG equations. § We will sketch the results of 
these techniques, referring the reader to [29, 30] for details. 

Unlike the excited state QPTs in nuclear physics described in Chap. 27, our 
QPT occurs in an inherently finite system. The QPT here disappears in the 
thermodynamic limit and cannot be extrapolated from finite size scaling or 
other such arguments. Thus there is no concept of nonanalyticity (see, e.g., 
Chap. 1) and it is formally a crossover. Given the prevalence of nano-devices 
and other inherently finite systems, we think the term finite-system QPT is 
a useful one, as we hope to convince the reader in the course of this chapter. 



5 An MPS approach is also possible, but requires spatial discretization, an aspect we wish 
to avoid. A continuum limit of MPS is possible in principle but tricky, requiring typically 
thousands of lattice sites; see [35]. 
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1.1 Fundamental Considerations 

We begin with a clear statement of the Hamiltonian. In position representa- 
tion the LLH is [4] 

#o = - £?=i d*. + g m E j<k m 6 k ), (1.1) 

where 0j is the azimuthal angle that satisfies < 6 3 < 2ir and giD is the binary 
contact interaction strength renormalized for our ID problem. We take the 
experimental context of the LLH as N bosons in a thin torus, or ring trap. 
The torus has radius R; we take the length, angular momentum, and energy 
to be measured in units of R, ft, and H 2 /(2mR 2 ), respectively, throughout our 
treatment. 

In a rotating frame of reference with an angular frequency 20, the LLH 
becomes 

H (O) = Ho - 2flL + Q 2 N , (1.2) 

where 

i' = -iE j =ide i (1-3) 

is the angular-momentum operator; we refer to Eq. (1.2) as the rotating LLH 
(rLLH). From the single- valuedness boundary condition of the many-body 
wave function [39], one can show that solving the eigenproblem in the rest 
frame i?o*o = E(Q)*f?o suffices in order to obtain solutions to the eigenproblem 
H(n)V = E(Q,)<£> [29]. The eigenvalue is then obtained from E(il) = E(0) - 
2Vt(L) + il 2 N, which is periodic with respect to fl. 

The subclass of solutions to the rLLH we focus on are the yrast states. 
This approach is useful for our ring system, because all the information about 
physical properties in the rotating frame are embedded within the spectrum in 
the rest frame. Thus the physical meanings of yrast states can be extracted by 
the simple transformation of yrast spectra. Since the LLH commutes with the 
angular-momentum operator, [Hq,L] = 0, the yrast problem is well defined 
irrespective of the sign and strength of interaction,^ and all the yrast states 
are eigenstates of both the LLH Hq and the rLLH H(fl). 

We can identify the yrast states by dividing the solution space of Ho and 
H into subspaces according to the two conserved quantities in the problem: 
the number of bosons N, and the total angular momentum L = (L). In each 
such subspace we can index excited states by q G {1, 2, . . .} in ascending order 
in energy; thus the state-ket is written \N,L;q). Then the yrast states are 
denoted as \N, L;q = 1). The essential properties of the ground and low-lying 



^Throughout this chapter we focus on repulsive interactions giD > 0, for which QPTs only 
occur in excited states. However, for attractive interactions, gm < 0, QPTs occur also in 
the ground state [31]. 
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excited states can be described within the yrast states, as we have verified 
explicitly with exact diagonalization studies." Thus we henceforth omit the 
quantum number q from the notations for eigensolutions. With the abbrevi- 
ation of the quantum number q = 1 and for fixed coupling constant <?id , the 
eigenvalues that correspond to the yrast states are written as En,l{^), where 
we explicitly write the parameter £1 in the notation in order to clarify in which 
frame the system is. With this notation, the eigenvalues in the nonrotating 
frame are written as £ , at i l(0). 

Most of the yrast states will turn out to be dark solitons for weak interac- 
tions, but in this limit a special subclass corresponds to uniform supcrflow. 
Due to the translational invariance of the LLH with respect to 9 and 0, prop- 
erties of a particular set of yrast states can be analyzed without solving the 
problem. These are the CMR states. In particular, they are the states for 
which the total angular momentum is equal to an integer multiple of the total 
number of atoms N. The energy of the CMR state takes the form 

E n , l =jn(0) = J 2 N + V int , (1.4) 

where Vint is the interaction energy and J 6 Z is an integer. We call J the 
center- of-mass quantum number, because it physically expresses the amount 
of uniform translation of the center-of-mass momentum. In the NLS mean- 
field theory, J is conventionally called the phase winding number, but this is 
a valid concept only in the weakly interacting limit, and we retain the more 
general notion for our purpose of identifying the QPT through all interaction 
regimes. In the rotating frame, the energy of the CMR state is given by 

E NJN (n) = (j-n) 2 N + v int , (i.5) 

where the change in energy associated with the frame change is involved only 
in the kinetic energy term, and the interaction energy is completely separated 
from the parameter $7. 

The ground state in the absence of the rotating drive is the state with zero 
angular momentum, -E/v,l=o(0)- The excitation energy of the CMR states 
with a finite angular momentum L = JN is thus given by 

E N ,j N (0) - E Nfi (0) = J 2 N, (1.6) 

which is independent of the strength of interaction gm . This is natural be- 
cause changing the total angular momentum by the amount JN leads to a 
Galilean transformation without changing the boundary condition for the or- 
der parameter. The ground state in the presence of the rotating drive is 
characterized by the CMR quantum number** 

Jo=L^ + l/2j. (1.7) 



II Such studies necessarily used a truncated single-particle angular-momentum basis, as the 

Hilbcrt space is otherwise too large. 

** \x\ is the floor function, meaning the largest integer that does not exceed x. 
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FIGURE 1.1 

Topological winding and unwinding of a BEC via formation of a dark soliton. 
Shown are the mean-field amplitude (solid curves with the left reference) and 
phase (dotted curves with the right reference). Uniform solutions with differ- 
ent values of the phase winding (i) J = 1 and (vi) J = are smoothly con- 
nected through the broken-symmetry dark soliton (ii)-(v) with a self-induced 
phase slip at (iv) $7 — 0.5. Reproduced from [28]. 



Because of the periodicity in the eigensolutions, an eigenstate \N,L) with 
the energy En.l (fi) has a denumerably infinite number of counterparts \N, L+ 
JN) and E N ^ L+JN (il), corresponding to arbitrary values of J € Z. Solving 
the yrast problem for a limited range of fixed angular- momentum states, e.g., 
—N/2 < L < N/2, therefore suffices to obtain all the eigensolutions. More- 
over, the spectra are degenerate for the same magnitude of angular momen- 
tum, -Ejv.l — -EV,-l in the absence of rotating drive, while this degeneracy 
is resolved in the presence of rotation due to the Sagnac effect [40] . All other 
yrast states for L out of this limited range can be obtained by shifting the total 
angular momentum by N while keeping the internal structure of the eigen- 
states. This is similar to a band theory concept [29], with —N/2 < L < N/2 
playing the role of the first Brillouin zone. 



1.2 Topological Winding and Unwinding: 
Theory 



Mean-Field 



We first introduce the basic concept of the finite-size metastable QPT in the 
simplest context, namely mean-field theory. We restrict ourselves to a fairly 
qualitative discussion; a complete quantitative theory, including full analytical 
solutions and perturbative stability analysis, can be found in [28, 29]. 
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The NLSE, which can be derived directly from the rLLH [1],^ is 

\{-id e - nf + g 1D N\m\ 2 ]m = HW), (1-8) 

where all bosons are taken to be in the same macroscopic mode which is 
described by the order parameter ip; this order parameter can be interpreted 
physically in terms of a density, |i/>| 2 7 and a phase, ip = Arg(ip). The single- 
valuedness of the wave function requires ip(6 + 2ir) = ip(0) + 2itJ, where J 
is an integer and can, in this weakly-interacting limit, be understood as a 
topological winding number. 

The uniform superflow solutions of the NLSE are just plane waves, ij>(0) = 
ipo exp(iJ#), while the dark-soliton train solutions are expressed in terms of 
Jacobi elliptic functions."'""'" We take the number of density notches in the dark- 
soliton train to be j. For the rest of this section, we consider the single soliton 
j = 1 for simplicity, but our discussion holds for arbitrary soliton trains J > 1 . 

Figure 1.1 illustrates how an initial uniform-superfiow solution with J = 1 
can be continuously unwound to J = 0. As tt increases starting from (i) the 
uniform superflow with J = 1, (ii) solitons start to form past a critical point 
^crit- ( ui ) The density notch deepens for Q^j> t < Q < 0.5. (iv) At Q = 0.5 
it forms a node, the phase of the soliton jumps by w, and the energies of the 
solitons with phase winding number 1 and are degenerate, (v) The soliton 
with phase winding J = deforms continuously as $1 increases, (vi) Finally, 
the state goes back to the uniform-superfiow state with phase winding J = 0. 
For the purposes of illustration we have taken a fixed interaction strength 
9idN = 0.6; in fact, interaction and/or rotation can be used to wind or 
unwind the order parameter; likewise, one can characterize the whole process 
in terms of system size, while holding other parameters fixed. 

In Fig 1.1 we observe the characteristic phase structure of dark solitons: 
inside the density notch, the phase changes more rapidly, and the total differ- 
ence across the notch yields the soliton's velocity; outside the density-notch 
region, the phase has a uniform slope, indicating the velocity of the uniform 
superflow. Since a dark soliton moves with a velocity which is opposite to 
the direction of increase in phase (to the left in Fig. 1.1), one can understand 
such stationary solutions as a cancelation between dark-soliton motion to the 
left, clockwise looking down on the ring from above, and uniform superflow 
to the right, or counter-clockwise. 

Figure 1.2(a) shows the energy difference between a dark soliton and uni- 
form superflow, as can be calculated analytically. This kind of bifurcation 
does not occur from the ground-state energy. However, for metastable states 
a bifurcation can occur between the uniform-superfiow state and the soli- 
ton state with the same winding number J. After bifurcation, the soliton 



ft The key assumption is to take a macroscopically occupied mode in the single particle 
density matrix and neglect small fluctuations around that mode; subsequently one replaces 
averages over products of operators with the product over averages. 
•f+In the infinite-system limit a single dark soliton is written in terms of a tanh function. 
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FIGURE 1.2 

(a) The bifurcation between unifbrm-superflow and dark-soliton states are 
illustrated in terms of (a) energy difference and (b) average angular momen- 
tum per particle L/N . The soliton solutions only exist in the area surrounded 
by two critical boundaries (white dotted curves) . The white arrow in (a) indi- 
cates the path taken during the unwinding of the phase shown in Fig. 1.1. In 
(b) one clearly sees two topologically distinct regions: in between the critical 
boundaries the characteristic quantization of L in the BEC breaks down, even 
though there is still macroscopic occupation of a single mode. Reproduced 
from [28]. 



energy becomes larger than the uniform-superflow energy. For convenience, 
Fig. 1.2(a) displays a white arrow indicating the higher-energy, soliton path 
taken in Fig. 1.1. Exact diagonalization studies were used to show that an 
arbitrary potential which breaks the symmetry of the ring automatically sets 
the system on this path. Shown in Fig. 1.2(b) is the key physical observ- 
able for our QPT, the average angular momentum per particle, expressed as 
L/N = J ddip* (-ide)ip, in the mean-field theory. Thus a continuous change of 
angular momentum is possible for ID Bose systems by taking the metastable 
states with energy slightly higher than that of the ground state. Both panels 
of Fig. 1.2 display a lower and an upper critical angular frequency, which we 



,(D 



>(2) 



call f^rit an d ^crif as a function of <?1D- 



Bifurcation of the soliton train from the uniform superflow constitutes a 
second-order QPT with respect to giD and/or ft. The derivatives of the 
soliton and uniform-superflow energies with respect to 51 have a kink at the 
boundary as can be verified analytically. This identifies the QPT, at least at 
the mean-field level, which occurs along a curve in the fl-gin plane. We also 
found that the Hessian for the energy is discontinuous along this curve, while 
the Hessian for the chemical potential diverges. 



10 Understanding Quantum Phase Transitions 



1.3 Finding the Critical Boundary: Bogoliubov Analysis 

We can better understand the QPT through the BdGE. The BdGE can be 
derived formally as the lowest order quantum fluctuations from the LLH; it 
can also be simply thought of as linear perturbation theory on the NLSE, and 
is taken as such in optics. In either case, without Bogoliubov analysis in terms 
of the BdGE we cannot be certain that our dark-soliton train solutions are 
stable even in the mean-field limit. Moreover, the BdGE allow us to explicitly 
identify the critical boundaries shown in Fig. 1.2. 

A stationary solution ip(9) of the NLSE under a small perturbation 5 evolves 
in time as 

${$, t) = e-^{V(0) + Y} 5u ^y~ lXnt + K(0)e iK % (1-9) 

n 

where (u n ,v n ), and A„ are the eigenstates and eigenvalues of the BdGE, 
respectively, and n £ |Z| denotes the energy-index of excitations. Recalling 
the structure of solutions in the Bogoliubov formalism [37], for each eigenvalue 
A n with positive norm, 

fS*d8[\u n (6)\*-\v n (6)\ 2 ]=h (1-10) 

there is also an eigenvalue A„ = — X n with negative norm. The BdGE predict 
an infinite set of such solutions. One exception to this rule can exist. This ex- 
ception corresponds to Nambu-Goldstone modes, sometimes called anomalous 
modes. For instance, consider a black soliton is at rest on the ring with respect 
to any background superflow in the rotating frame. The Goldstone mode of 
the soliton corresponds to the soliton moving at a small constant velocity in 
this frame. We need consider only the Goldstone mode and eigenstates that 
satisfy (1.10), since the eigenvalues with negative norm do not have physical 
meaning. For the Goldstone mode the corresponding eigenvalue is zero, and 
the latter eigenstates that have positive norms can be real (i.e., positive or 
negative) or complex depending on the stability of the condensate mode. In 
general, A„ <G C. 

The excitations of a uniform superflow with phase-winding J are straight- 
forward to calculate [37]. Fluctuations from that condensate mode are given 
by the eigensolutions of the BdGE with positive norm, 

X\ J ^ = y /p(p + 2^)-2l{n-J), «,(xe ! ( J + ! » 9 , vi oc e - i( - J -W, (1.11) 

where the 'us' subscript stands for uniform superfluid, 7 = giYyN/2ir, and 
I G Z denotes the single-particle angular momenta of the excitation, a good 
quantum number since [H(Cl),L] = 0. 

For 7 > 0, all the eigenvalues A ; ' are real, as apparent from Eq. (1.11). 
From Eq. (1.11) we also find that several negative eigenvalues, associated 
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with eigenstates of positive norm, appear when we take a metastable excited 
state as a condensate mode. These negative eigenvalues correspond to other 
plane-wave branches located in lower energy regimes than the input conden- 
sate mode itself. For the case of repulsive interactions, the number of negative 
eigenvalues thus coincides with the number of stationary states that are lo- 
cated in a lower energy regime than the metastable state under consideration. 
Solutions of the BdGE for dark-soliton trains require a numerical calculation. 
We can identify the critical boundary (7 cr it, ^crit) for the QPT by finding 
where a particular negative eigenvalue changes its sign. One equates I to be 
±j in Eq. (1.11) and imposes the condition on the eigenvalue, Aj_'^ = 0. 
Then the critical boundaries are given implicitly by 



ficrit" J = ±V / 072) 2 +7crit/2, (1.12) 

where the ± sign gives the lower and upper boundary, respectively, which we 

called £l„lt ana - ^crit- The region 7 < 7 cr ;t is identical to requiring that the 
critical-boundary condition has a real solution. 

In Fig. 1.3 we illustrate the QPT from uniform-superflow to dark-soliton 
solutions; we illustrate only < < 0.5 since the eigenvalues are symmetric 
for the other half of the Brillouin zone. The soliton-train solutions occurring 
for fl > f^ rit are therefore linearly stable.* The excitation energies from the 
soliton branch are found to be close to those from the plane- wave branch. The 
notable feature in the soliton regime is that there appears a Nambu-Goldstone 
mode, which is continuously connected with one of the negative eigenstates 
with I — Sj from the plane wave. This mode reflects the spontaneous symme- 
try breaking of the soliton-train state. At the point f2 no dcs — 0.5, a degenerate 
pair of excitation branches emerges, where the phase jumps up or down by ir 
at each soliton in the soliton train. This degeneracy is similar to that seen for 
a single soliton in Fig. 1.1, but slightly complicated by the multiple solitons 
and higher phase winding. For £1 > f2 n odes the excitation branches from the 
soliton train are also symmetric with respect to f2 no des- 

For expediency we skip more details of BdGE solutions and summarize 
the overall results. We found that excitations from the uniform supcrflow in 
the j th excited-state have j thermodynamically unstable modes within the 
critical region in the (7, fi)-plane. In fi„ it < fl < 0.5 in the dark-soliton train 
states become stable. At the critical boundary, all energies of the stationary 
solutions and the eigenvalues of the BdG equations continuously connect to 
the excitations in the soliton regime without any energy discontinuity. In 
particular, when one of the excitation energies from a plane-wave metastable 
state changes sign from negative to positive, a soliton branch with the same 
phase- winding number appears, and the Nambu-Goldstone mode manifests as 
a result. 



*They are also nonlinear ly stable [33]. 
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FIGURE 1.3 

Bogoliubov analysis of uniform superfiow (thin dashed curves) and dark- 
soliton trains (thick solid curves). The initial uniform-superfiow state (zero 
excitation energy shown as the thin solid curve and labeled with text) with 
phase winding J = 2 becomes the Nambu-Goldstone mode of the three-notch 
(j = 3) dark-soliton train, past the critical frequency ^rit- Other curves show 
excitation energies from the uniform and soliton train states, respectively, with 
various angular momenta labels /. We take fixed interaction gn)N /2ii = 1 for 
this illustration. Reproduced from [29]. 



Thus, upon consideration of general phase winding J and soliton-number 
j, one realizes that there are a denumerably infinite set of paths to connect 
uniform-superfiow states via soliton trains in a metastable system of scalar 
bosons on a ring. Associated with this transition, the energy of the solitons bi- 
furcates, and a continuous change in the angular momentum becomes possible 
in the mean-field theory past a critical boundary in the interaction-rotation 
plane. In general, bifurcations can indicate symmetry breaking associated 
with a QPT, as in the elementary example of the quadratic-quartic potential, 
where a single potential well bifurcates into two as the relative strengths of 
quadratic and quartic terms are tuned. 
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FIGURE 1.4 

Yrast energy eigenstates for N — 40 bosons on the ring, obtained by exact 
diagonalization of the rLLH with the cutoff angular momentum |Z C | = 2. 
The spectrum has a kink where L is an integer multiple of N . Everything 
is symmetric for L < 0, and in the large N limit, the number of points 
increases and the discrete yrast energies approach a continuous curve while 
the curvature and kink points remain unchanged. Reproduced from [30]. 



1.4 Weakly-Interacting Many-Body Theory: Exact Di- 
agonalization 

We now show how in general to distill the uniform-superflow and soliton- 
train mean-field branches from a sea of many-body eigenvalues. We show 
how the mean-field soliton branch, for which average angular momentum is 
not quantized as an integer, emerges from the yrast spectra. The meaning 
of spectra related to symmetry breaking and the Nambu-Goldstone mode 
associated with the existence of the soliton branch is also discussed. 

In order to take the next step towards a microscopic many-body solution 
of the LLH, it is convenient to rewrite it in second-quantized form. One 
expands the bosonic field operator in terms of a plane-wave basis with the 
single-particle angular momentum Z, 

m^i^y^ZtZ^e* 16 , (1.13) 

where the pre-factor of (2tt)~ 1 ' 2 comes from the normalization of the plane 
wave, and bi and bj are annihilation and creation operators which obey the 
usual commutation relations for bosons. Equation (1.13) satisfies the periodic 
boundary condition ip(9) = ip(0 + 2ir). Then one finds 

^0 = Et^oo l 2 bfh + fflD E2Tm,n=-oo &$6mU*+J ltn+ n • (1-14) 

The eigenstates can be expanded in terms of a Fock-state basis |{nz}) that 
represents the occupation number of each single-particle angular-momentum 
state, 

\{ni}} = |...,n_i,n ,ni,...). (1.15) 
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FIGURE 1.5 

(a) Performing a Legendre transform on the yrast states in Fig. 1.4, we 
obtain the full yrast spectrum of the rLLH as a function of fl. Each curve is 
distinguished by a different total angular momentum L. (b) A zoom of (a) 
in the important region where solitons appear in the mean-field theory. The 
dashed curves show a comparison to the NLSE superflow (blue) and soliton 
(red) branches. The enclosed region is called a swallowtail. Reproduced 
from [301. 



These states satisfy the conservation laws 



(1.16) 



For numerical calculations we use a cutoff angular momentum l c > 0; thus 
L e [—l c N,l c N]. In the weakly interacting regime gmN ^ <^(1)> a cutoff of 
l c = 2 provides a quantitative agreement in energy eigenvalues^ with those 
obtained by the Bethe ansatz shown in Sec. 1.6. 

Figure 1.4 shows the yrast energies for the nonrotating LLH, En.l(£1 = 
0) = (N,L\H \N,L), for interaction strength g 1D = 2.5 x 10~ 2 7r and number 
of bosons N — 40. The ratio of the mean-field interaction energy to the kinetic 
energy corresponding to these values of giD and N is gir,N/(2n) = 0.5. The 
kink observed in the figure is exactly the CMR state L = JN, here L = N. 
In the mean field, this would be the uniform-superflow solution. The other 
non-CMR yrast states contribute to the dark-soliton solutions, as we will 
show. 

To generalize Fig. 1.4 to the rotating rLLH, all we have to do is follow the 
procedure laid out in Sec. 1.1. We perform a Legendre transform on the LLH 



tThe convergence is better for (i) a larger number of atoms, (ii) lower eigenvalues, and (iii) 
smaller strength of interaction. Figure 1.9 indicates that for N = 10 and up to the yrast 
states \L\ = 2N, the deviation of the exact diagonalization results from the Bethe ansatz 
results starts to emerge for gio > 1. The results for N = 40 in this section has, in general, 
better agreement with the rigorous results. 
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yrast spectrum, 

E N . L (0) -4 E NiL (il) = E NtL (0) - 2Q.L + n 2 N. (1.17) 

Figure 1.5 plots transformed yrast energies En.l(Q), where the finite number 
of points from Fig. 1.4, each of which was characterized by a different angular 
momentum L, become convex downward curves in Fig. 1.5. Each curve is thus 
characterized by a different total angular momentum and has a minimum at 
a certain value of fl. 

The energy En^l = j n(Q), shown in Fig. 1.5, corresponds to the ground 
state where Jo is the ground-state CMR quantum number given by Eq. (1.7). 
The angular-momentum states with L — JN correspond to the CMR states, 
and the center of the parabola is located at J7 <G {Z} at which the CMR state 
becomes the ground state. 

In Fig. 1.5(a) we observe an extremely high density of states around O € 
{±0.5, ±1.5, . . . } due to the crossing of many eigenvalues, as shown more 
clearly in the zoom in panel (b). The region of high density takes the shape 
of a swallowtail; the same shape was found to occur purely within mean-held 
theory, past the critical boundary for the QPT [29]. This swallowtail region 
is almost hlled by various energy eigenvalues of various angular-momentum 
states crossing each other. The domain with the high-density swallowtail 
shape looks as if it is enclosed by the two kinds of stationary branches pre- 
dicted by the mean-held theory. 

How can we understand this spectrum? It is an interesting fact that mean- 
held theories tend to break the symmetries of many-body theories. In our case, 
if the NLSE has a soliton solution, the solution must be localized on the ring; 
the NLSE is nonlinear and cannot work with superpositions of solitons. In 
contrast, the quantum theory hnds a superposition of all possible positions. * 
This is the main reason why Dziamarga et al. found that an initial mean-held 
soliton delocalizes at a higher level of quantum theory [14]. 

We can also see this effect in vortex formation in a scalar condensate under 
rotation. Solving the yrast problem in 2D results in all the angular-momentum 
states, including the at- rest condensate (L = 0), off-axis vortex (0 < L < 
N), a centered vortex (L = N), and vortex lattices (L > N). However, 
in experiments one drives the system with a specihe angular frequency. In 
such a situation, there exists a small distortion in the trap, which selects 
a metastable angular-momentum state with respect to the variation in the 
angular momentum of the condensate. As a result, in reality one does not 
observe a stationary off-centered vortex except as a transient state. 

The same argument applies to our case. In the presence of any kind of 
noise, such as an infinitesimal distortion of the trapping potential, quantum 
measurement of the matter wave, or whatever else breaks the translation sym- 
metry of the ring trap, the realizable stationary state or metastable stationary 



•f There is a quantum entropy associated with this degeneracy [32], 
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FIGURE 1.6 

The mean-field solutions can be extracted from our microscopic many-body 
analysis by using an extremization condition, (a) Energy. The thin curves 
have integer average angular momenta and agree with the energy of the mean- 
field uniform-superflow states; the thick points signify nonintcgcr angular mo- 
mentum and have an energy close to that of the mean-field dark-soliton states, 
(b) Angular momentum. Note that we use the same parameter set as Figs. 1.4 
and 1.5. Reproduced from [30]. 



state is determined by extremization with respect to variations in angular mo- 
mentum. In order to find the metastable states we impose the condition 



d L E NjL (n) = o. 



(1.18) 



with $7 and <?id held fixed. Figure 1.6(a) plots energy eigenvalues that satisfy 
this condition as a function of $7; and Fig. 1.6(b) shows the corresponding 
angular momentum. These figures are quite similar to those given by mean- 
field theory, i.e., by imposing the stationary condition (1.18) for the manifold 
of eigenvalues we identify the mean-field stationary branches. 

Referring to Fig. 1.6(a), we cannot call the thick curve a soliton branch in 
a rigorous sense, because each point is an eigenvalue of the rLLH and thus 
the associated eigenstate still possesses translational symmetry, unlike mean- 
field dark solitons. Instead, we call all the angular-momentum states inside 
the swallowtail in Fig. 1.5 the soliton components, because in the presence 
of infinitesimal noise these states do form a broken-symmetry state, which 
we denote |%). Soliton solutions of the NLSE can be interpreted in terms of 
the eigensolutions of the rLLH as a state where the several eigenvalues in the 
swallowtail region are collectively superimposed. 

The energy associated with this superposition does not change significantly 
because the energy required to make it is on the order of 1/JV. As a result, the 
energy of the broken-symmetry soliton state |%) is also well approximated by 
the thick curve in Fig. 1.6(a). In the presence of an infinitesimal symmetry- 
breaking potential, the angular momentum is no longer a good quantum num- 
ber. However, the expectation value of the angular momentum (x|£|x) agrees 
well with that of the dark solitons obtained by mean-field theory, and thus 
behaves like that shown in Fig. 1.6(b) [28]. 
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With all these caveats in mind, we briefly state that the branch drawn by 
the thick curve in Fig. 1.6(a) is the quantum soliton branch in the weakly 
interacting regime. 

We also calculated the second derivative d 2 E^ i L(n)/dL 2 with respect to fl 
in order to check whether the metastable angular-momentum state is a local 
maximum or minimum. The superflow state with a CM quantum number 
Jo = [ft + 1/2J is indeed the ground state because the second derivative is 
positive at that point, while the thick points are local maxima with respect 
to L, since the second derivative is negative. 

Let us summarize our first foray into the many-body theory. In this section 
we obtained the yrast states \N,L) of the rLLH by exact diagonalization. 
We used an extremization condition to extract the mean-field results, among 
them the characteristic swallowtail shape in the spectrum. We again found 
two kinds of metastable branches, uniform superflow and quantum soliton, 
consistent with mean-field theory. The region of a high density of states, 
where the different angular-momentum states in the quantum theory densely 
cross, indeed agrees with the soliton regime predicted by the mean-field and 
Bogoliubov theories. The phrase "quantum flesh sewn onto classical bones" 
has been used elsewhere [41] as a visual metaphor, perhaps inspired by x-ray 
images, to describe this accord. 



1.5 Strongly-Interacting Many-Body Theory: Tonks-Girardeau 
Limit 

Let us now move to the other extreme from the mean field. Taking the 
interactions to be very large, the bosons become impenetrable and thus can 
be mapped to spinless fermions. This is the TG regime. We again solve the 
yrast eigenproblcm and study the consequences. In particular, we introduce 
the particle and hole excitations, which are well defined in the fermionized 
gas, and show that these excitations are related to the mean-field stationary 
states in the opposite weakly interacting limit. 

To be precise, the TG Bose-Fermi mapping theorem [38] is 

* B ({e}) = l[sgn(d i -d j )* F ({e}), (1.19) 

i>j 

where ^b({^}) is the bosonic many-body wavefunction, expressed in first 
quantization, and typ({9}) is the equivalent spinless fermionic many-body 
wavefunction. This theorem holds for all the eigensolutions [42], and hence 
significantly simplifies our eigenproblcm.^ 



' For simplicity of notation we will show the analytic expression only for an odd total number 
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FIGURE 1.7 

(Left) A sketch of how to construct TG yrast states. (I) Particle excitations 
where the angular momentum of a particle increases while a hole is positioned 
at lp + J. (II) Hole excitations where a particle is placed at the lowest 
unoccupied state and the angular momentum of a hole decreases. When L is 
an integer multiple of N, only the center-of-mass angular momentum is shifted 
from the ground-state configuration. (Right) Type I and type II excitation 
energies as a function of L/N for N — 11 free fermions. Reproduced from [30]. 



The ground state of N (odd) free fermions is obtained by the occupation of 
the angular-momentum states from I = —lp to / = lp [sec L = in Fig. 1.7], 
where lp = (N — l)/2 is the Fermi momentum. The ground-state energy is 
thus 



E N , L=0 (n = 0) = Ei=-i F l 2 = N(N 2 - 1)/12. 



(1.20) 



Beyond the ground state, Lieb has shown [3] that excitation of the repulsively 
interacting Bose gas in the thermodynamic limit has two branches. The first 
branch is called type I and was shown to be in agreement with the Bogoliubov 
spectrum of plane waves in the weakly-interacting regime. The second branch 
is called type II, and this was supposed to be absent in the Bogoliubov spec- 
trum. For the fermionic formulation of the TG gas we observe that the type 
I and II branches correspond to particle and hole excitations, respectively. 

Let us reconsider these branches in the context of yrast states. For the 
excited state En^(0) we can find the particle and hole excitations. The 
procedure is slightly technical, but we provide a sketch in Fig. 1.7(a), and a 
brief description below. 



of particles. For an even number of particles the periodic boundary condition must be taken 
as antisymmetric. 
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Type T. Remove a particle at the Fermi momentum lp and place it at the 
momentum lp + L. For free fermions, there is no energy-level reconstruction 
in an (JV± l)-particle system associated with removal or addition of a particle. 
The energy of the type I excited state E N L (0) is thus obtained as E N L (0) = 
^jv,o(0) — lp + {If + L) 2 , or, relative to the ground state, 

A£$ L (0) ee E$ iL (0) - E Nfi (0) = L(N + L- 1). (1.21) 

There is no limitation on the single-particle angular momentum for Type 

I excitations. Such excitations are doubly degenerate for O = 0, for I — > 
—I. The excitations from the CMR states are similarly obtained: to get the 
excitations from the J = JN state, remove a particle at lp + J and replace 
it at I — lp + J + L — JN. The resulting excitation energy is given by 
AE$] L (0) = J 2 N + (l F + J + L - JN) 2 - (l F + J) 2 , L > JN. 

Type II: Starting from the ground state, remove a particle (create a hole) at 
the momentum lp — L+1 and place the particle at lp + 1, where < L < N. It 
is clear from Fig. 1.7 that the hole with this kind of low- lying excitation energy 
be created only within the range — lp < I < lp. The energy of this excited 
state is given by £^^(0) = E N , L (0) - (lp - L + l) 2 + (lp + 1) 2 , or, relative to 

the ground state, AE^ L (0) = L(N -L + l),0<L<N. Starting from the 
CMR state L = JN, we can extend this procedure to L > N, producing the 
series of humps shown in Fig. 1.7. From the fact that the excitation energy 
of the CMR state L = JN is AEn^ = jn(0) = J 2 N, we can obtain the type 

II hole excitation energy for the general case as 

AE^ L (0) = J 2 N +{l F + J+l) 2 -{l F -L + JN + J+ l) 2 , (1.22) 

where JN < L < (J + 1)N and AE^> L (0) = E^ L (0) - E Nfi {0) is again 
the energy relative to the ground state. We note that Lieb's original study 
focused on the region < L/N < 0.5 [3]. 

Next we rotate these states to the TG yrast states for the rLLH, again 
using our Legendre transform. The energy of the type II excited state for 
JN < L < (J + 1)N measured relative to -Eat„l=o(0) is given by 

£jjg.(fi) = AE { ^ L (0)-2nL + n 2 N. (1.23) 

We can again get the key solutions from the extremization condition (1.18), 
where now we take the energy to be that of type II excitations. By inspec- 
tion, CMR states L = JN are metastable states, either ground or excited. 
Condition (1.18) gives another metastable angular momentum, 

L = (N + l)(J+l/2)-n, (1.24) 

and the corresponding energy, 

EflQl) = (N + l)[fi - (J + 1/2)] 2 + N(N + l)/4. (1.25) 
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FIGURE 1.8 

We can again use an extremization condition to get the realizable yrast 
eigenstates in the strongly-interacting TG limit. Shown is the result obtained 
from the TG Bose- Fermi mapping for N — 11 bosons. All energies are plotted 
relative to the nonrotated ground-state energy. The thin curves show the 
neighboring CMR states with angular momentum L — JN and L = ( J+ 1)N. 
The bold curve shows the intervening type II states with angular momenta 
between these values. The angular momentum decreases from (J + 1)N to 
JN according to Eq. (1.24) along the bold curve. The arrows show the critical 
frequencies for the QPT. Reproduced from [30]. 



As a function of fi this is a parabolic curve, shown in Fig. 1.8. 

We again find certain critical angular frequencies where the metastable type 
II branch disappears and merges into the CMR branch: 



^ r = (J+l/2) T iV/2. 



(1.26) 



The stable angular momentum approaches L = ( J + 1) N at £l~ and L = JN 
at fi^., respectively, and the corresponding energy coincides with the energy 
of CMR states. Equation (1.26) determines the critical boundary for our QPT 
in the strongly interacting limit. 

This whole picture closely matches that of the uniform-superflow to dark- 
soliton transition in the weakly interacting limit, where there exists a critical 
angular frequency at which the soliton branch bifurcates from the superflow 
branch: here type II yrast states branch off of the type I CMR states. We 
therefore interpret the hole excitations in the TG limit with the soliton branch 
in the weakly interacting limit. It remains to fill in the picture between, as 
we proceed to do in the following section via the Bethe equations. 
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FIGURE 1.9 

Evaluation of three many-body solution methods, Bethe equations, exact 
diagonalization, and mean-field NLSE (labeled GP Theory), in the weakly- 
to medium-interacting regime in terms of yrast eigenenergies (note log scale) . 
We take N = 10 particles and CMR quantum number J G {0, 1, 2} only for 
the purposes of illustration. The inset enlarges the medium-interacting regime 
on a linear scale. Reproduced from [30]. 



1.6 Bridging All Regimes: Finite-Size Bethe Ansatz 

In this section we will show that our QPT stretches over the whole interaction 
range, and also validate our physical interpretation of Lieb Type II excitations. 
We note that, complementary to our work, the spectrum of the LLH has been 
obtained [43] by treating the inverse of the TG parameter, which is infinite at 
the TG regime, as the expansion parameter, and its analytical interpolation 
was given recently [44]. 

Using the Bethe ansatz, we obtain the N simultaneous Bethe equations for 
our system, 

\iV-l -2irie n _ TT *n ~ *m + *31d/2 



(-1) J 



n 



£ n 



igm/2' 



(1.27) 



which determine the set of values {^ n }, called quasi-angular momenta for 
each atom n G {1, ..., JV}." The quasi-angular momenta fully characterize 
the angular momentum and energy, given by L — J2n=i^n an< ^ En,l(Q = 
°) = En=i4> respectively. 

We numerically solve the real part of the Bethe equations (1-27) for each 
set of energy levels characterized by the different total angular momenta. The 



^It is sufficient to solve for only the real part of these equations [30]. 
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FIGURE 1.10 

(a) Solutions of the Bethe equations in the nonrotating frame for a range of 
angular momenta \L\ < 2N and N = 10 particles. The horizontal dashed line 
corresponds to the energy of free fermions (1-20) with the same number of 
atoms, (b) Excitation energies of yrast states as a function of average angular 
momentum \L\/N for fixed strengths of interaction. A, B, C, D, and E are 
marked in (a). The kinks are again the CMR states. Reproduced from [30]. 



numerical solution of Eqs. (1.27) is highly sensitive to the initial set of trial 
values of {in}- If this initial set is sufficiently close to a solution for a target 
angular-momentum state, the set of solutions {£ n } can be correctly obtained. 
In contrast, if the initial set is closer to another angular-momentum state, the 
total angular momentum reveals undesired jumps, deviating from the target 
angular momentum. In such a case we again start from another initial set of 
trial values of quasi- moment a. Our algorithm is outlined in [30] . 

In Fig 1.9, we first perform a comparative study for mean-field, TG, and 
Bethe solution methods in the weakly- to medium- interacting regime." As 
can be seen, for #id < 0(1), the three methods agree to a few percent. 

In Fig. 1.10(a) we show the spectra over the entire interaction range, for 
the system in the nonrotating frame: specifically, En,l(Q = 0) for \L\ g 
{0, ±1, . . . , ±2N}. All ground- and excited-state energies monotonically in- 
crease with respect to giu > 0. However, the energy does not monotonically 
increase with respect to the total angular momentum for a fixed strength 
of interaction. This is illustrated in Fig. 1.10(b). The labels A, B, C, D, 
and E correspond to the cross-cuts labeled in panel (a). While in the very 
weakly-interacting limit (curve A) the spectrum still looks almost linear, as 
the interaction increases, the kinks in the yrast spectra at the location of 
CMR state, L — JN, become more pronounced due to the large increase in 



1 1 We note that the concept of yrast state for the angular momenta JN < L < (J + l)iV does 
not exist in the mean-field theory: this theory is concerned only with the single-particle 
angular momentum, which coincides with the average angular momentum in this theory. 
We thus plot the mean-field energy for the integral single-particle angular momenta. 
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FIGURE 1.11 

Spectra of the rLLH obtained via the Bethe equations and satisfying the ex- 
trcmization condition for (a) weakly-interacting regime, (b) and (c) medium- 
interacting regime, and (d) strong-interacting regime. The thin curves show 
CMR states, while the thick curves show yrast states with L/N noninteger, 
and come from the type II excitation branch. Reproduced from [30]. 



the energy of the yrast states with L in between neighboring CMR states 
at JN and (J + 1)N. For strong interactions (curve E), the system is in 
the TG regime; compare to the TG calculation for Type II excitations from 
Fig. 1.7(b). We observe numerically that the excitation energy of the CMR 
state L = JN is independent of girt, and is given by Eq. (1.6), namely 
{E n ,l - E Nfi )/N = (L/N) 2 . This follows from the nature of the CMR state 
L = JN, which is just a Galilean boost of the nonrotating state; under this 
transformation interactions are unchanged. 

We proceed to the rotating frame, as before, via our Legendre transform, 
which can be framed in terms of quasi-angular momenta as 



EnM = Ef= 



^) 2 = E N}L (0) - 2flL + Q 2 N. 



(1.28) 



The results are shown in Fig. 1.11 for various strengths of interaction. The 
thin curves are parabolas (O — J) 2 + Vint for various values of center-of-mass 
quantum numbers J. The lowest possible energy of the CMR state is thus 
given by Vint at 17 = J G {Z}. The thick curves plot other stable yrast states 
of noninteger L/N from the condition (1.18). The weakly-interacting mean- 
field regime is shown in Fig. 1.11(a), where the type II branch that satisfies 
the metastable condition just starts to appear. Thus these are the energies 
of the quantum solitons [see also Fig. 1.6(a)]. As the interaction increases 
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FIGURE 1.12 

(a) Two CMR states with L = 0, N (thin curves) and the type II branch (thick 
curve) that connects them. The energy is defined relative to the interaction 
energy V; nt . (b) The quantum phase diagram over all interaction strengths; in 
the gray region yrast states corresponding to Lieb's type II excitations allow 
for a continuous change in angular momentum, while in the white region only 
CMR states with quantized L/N are allowed, (c) Enlargement of (b) in the 
weakly-interacting regime. The solid curve is the critical boundary given by 
the Bogoliubov theory, (d) Difference between the phase boundaries given 
by the Bethe ansatz and the Bogoliubov theory. The Bogoliubov critical 
boundary overestimates the correct Bethe result. Reproduced from [30]. 



[Figs. 1.11(b) and 1.11(c)] the domain with the swallowtail shape enclosed by 
the two CMR branches, as well as the size of the type II branch, increases. In 
the TG limit [1.11(d)], the area of the swallowtail region saturates the spectra. 

These behaviors are quantitatively summarized in Fig. 1.12(a), which shows 
the energy En,l{Q)/N of metastable states relative to the interaction energy 
Vint at each strength of interaction. The CMR branches drawn by the thin 
curves no longer have a c/iD-dependence because of the subtraction of V5 n t, 
while the thick curve gradually increases the domain over which it extends as 
the interaction increases. For simplicity we plot only two CMR branches with 
angular momenta L — 0, N, and a metastable state associated with the type 
II branch that smoothly connects these two CMR states. 

As fi increases, the thick curve bifurcates from the CMR branch with an- 
gular momentum L = N at a certain critical f2 < 0.5, and at £1 — 0.5 the 
energy is minimal. As Q increases further, this branch smoothly merges into 
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the CMR branch with angular momentum L — and eventually disappears 
at a certain critical ft > 0.5. We therefore find that the same kind of energy 
bifurcation which was found in the mean-field theory persists over the full 
range of repulsive interactions. Figure 1.12(b) illustrates this idea by plotting 
the existence range of the metastable-state type II excitation branch. The 
shaded area indicates the existence of such a branch. The angular momentum 
on the lower critical boundary is given by L — TV, and the angular momentum 
linearly decreases as ft increases, just like in Fig. 1.6(b). At O = 0.5 (ver- 
tical dashed line), the value of the angular momentum is given by L = N/2 
irrespective of the strength of interaction. At a certain value of il(> 0.5) the 
angular momentum eventually goes to zero, causing the metastable hole exci- 
tation branch to disappear. This behavior corresponds to the fact that in the 
mean-field theory the type II branch bifurcates from the uniform-supcrflow 
regime, developing nodes, and it again merges into the uniform-supcrflow 
regime with the increase of Q [28, 29]. The critical boundary approaches 
51 = ( J + 1/2) ± N/2 in the strongly-interacting regime. 

Finally, we explicitly evaluate the BdGE analytical predictions for the criti- 
cal boundary. In Fig. 1.12(c)-(d) the Bethe-equation prediction for the critical 
boundary is compared with that obtained from the BdGE in the weakly- 
interacting regime. We observe that Bogoliubov theory predicts the quanti- 
tatively correct critical boundary to the 5% level up to gm < 5 (for N = 10), 
but it significantly overestimates the boundary as the interaction increases. 

Thus we find a critical boundary for all values of interaction strength. 



1.7 Conclusions and Outlook 

We have tied together a number of outstanding issues in ultracold quantum 
gases. The first is the understanding of quantum effects on characteristic 
mean-field solutions, namely dark solitons. The second is the concept of a 
phase transition in a mesoscopic system which is strongly isolated from its 
environment. Other issues we addressed included the physical interpretation 
of Lieb's Type-II excitations and the correspondence of various ID techniques, 
ranging from the NLSE to the Bethe equations to the TG mapping. 

The first issue has occupied many people over the last decade. Work in this 
area includes quantum derealization of dark solitons in the weakly interact- 
ing regime [14, 19], decay of dark solitons in the more strongly interacting 
regime [15], the search for dark solitons in the extreme limit of the TG gas 
via the TG Bose- Fermi mapping [46] , and most recently the role of thermal 
fluctuations in defect formation by the Kibblc-Zurck mechanism [21]. Our 
results indicate that while the mean-field dark soliton bears some meaning in 
the weakly-interacting regime, where one can make a close correspondence be- 
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tween the underlying many-body quantum soliton branch and the symmetry- 
broken mean-field soliton, in the medium- to strongly-interacting regime no 
such correspondence is likely. Thus for stronger interactions we do not expect 
to observe a density notch with a characteristic phase structure. Given recent 
experiments tuning the interactions over seven orders of magnitude and the 
increasing commonness of ring traps, we hope to sec tests of these ideas in the 
near future. A new dynamical study by Brand and Kolovsky [45] indicates 
that in various regimes stretching from the NLSE to the TG gas our ideas 
have some merit; after all, if our predicted many-body generalization of dark 
solitons took the lifetime of the universe to nucleate, our ideas would remain 
restricted to the realm of theory; this appears not to be the case. A comple- 
mentary study to Brand and Kolovsky's might use MPS methods to address 
this question. 

The second issue is trickier. We all have a thermodynamic view of phase 
transitions and states of matter drummed into us, and new contexts such as 
graphcnc arc showing that some of these ideas have to be generalized. In our 
case we find a QPT, technically a crossover, that nevertheless displays char- 
acteristics of phase transitions: topological properties of the system change 
and the energy reveals a cusp at some level of derivative. Moreover, our QPT 
is inherently finite size and has no meaning in the thermodynamic limit. For 
isolated systems such as ultracold quantum gases or nuclei (see Chap. 27) a 
QPT in excited states is a useful concept, in contrast to QPTs in solid-state 
materials. The mean-field theory is expressed in a nonlinear partial differen- 
tial equation, the NLSE; it is useful to note that bifurcations appearing in the 
solution space of such effective nonlinear theories often point to a QPT in the 
underlying linear many-body theory. 

As a very brief summary of other issues addressed in this chapter, we showed 
that the soliton-train uniform-superflow picture falls naturally out of consider- 
ation of center-of-mass-rotation states and quantum soliton states, the latter 
being yrast states contained within a swallowtail. These considerations pro- 
vided a complete description of the finite-system metastablc QPT from the 
weakly- to strongly-interacting regimes. In the process we cleared up the 
long-standing problem of the physical interpretation of Lieb's type-II or hole 
excitation branch: type-II excitations are exactly the yrast states, which in 
the weakly- interacting limit can be expressed as dark-soliton trains, appearing 
in broken-symmetry form in mean-field theory. 

Our static study can be extended to many different contexts in ultracold 
quantum gases, for example, spinor Bose systems, where the hyperfine struc- 
ture of the constituent atoms in the BEC plays the role of a spin, and one 
finds spin-one or even spin-two and higher models, as described in Chap. 10. 
A number of metastable and/or finite-size QPTs are possible in these systems. 
Ultracold fcrmions also provide a fine candidate, and we expect a metastable 
QPT occurs in the Bardeen-Cooper-Schricffcr (BCS) superconducting phase. 
One could trace this QPT through the BCS-to-BEC transition. This prob- 
lem is particularly interesting halfway through the crossover in the unitary 



Metastable Quantum Phase Transitions in a One- Dimensional Bose Gas 27 

regime,** a many-body system under intensive investigation in ultracold quan- 
tum gases, in the quark gluon plasma (see Chap. 25), and in the AdS/CFT 
(anti-de Sitter / conformal field theory) mapping (see Chap. 28). 

More generally, we can ask if the concept of QPTs in finite systems is useful 
beyond ultracold quantum gases. What happens to critical exponents? Can 
they be defined without the use of finite-size scaling? Is the concept of a 
universality class useful? Can backing off from the thermodynamic limit but 
retaining the QPT concept in some fashion help us better understand the 
fundamental physics of mesoscopic systems and help us invent better nano- 
devices? 
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